Incoherence of Bose-Einstein condensates at supersonic speeds due to quantum noise 
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We calculate the effect of quantum noise in supersonic transport of Bose-Einstein condensates. 
When an obstacle obstructs the flow of atoms, quantum fluctuations cause atoms to be scattered 
incoherently into random directions. This suppresses the propagation of Cherenkov radiation, creat- 
ing quantum turbulence and a crescent of incoherent atoms around the obstacle. We observe similar 
dynamics if the BEC is stirred by a laser beam: crescents of incoherent atoms are emitted from the 
laser's turning-points. Finally, we investigate supersonic flow through a disordered potential, and 
flnd that the quantum fluctuations generate an accumulation of incoherent atoms as the condensate 
enters the disorder. 

PACS numbers: 



CD 



I 

o 
o 



> 

00 

d 

00 

o 

> 

X 
J3 



I. INTRODUCTION 

The production of elementary excitations above a crit- 
ical velocity is an ubiquitous feature of fluids. A fa- 
miliar example is the sonic boom created as an air- 
craft exceeds the speed of sound in air. Recent exper- 
iments py have demonstrated the analogous effect in di- 
lute gas Bose-Einstein condensates (BECs): Cherenkov 
radiation which precedes an obstacle travelling at super- 
sonic speeds relative to the BEC. This phenomenon has 
been reproduced theoretically 0, Q by numerical inte- 
gration of the Gross-Pitaevskii (GP) equation. However, 
unlike the aircraft, the BEC is a quantum object, so we 
might expect quantum fluctuations to play a significant 
role in the dynamics, especially because quantum noise 
has been shown to generate scattering halos and quantum 
turbulence for large collisional velocities Moreover, 
superfluids differ from normal fluids in that they have 
a critical velocity, known as the Landau critical veloc- 
ity , for the onset of dissipation. One of the early tests 
of superfluidity in BECs was to demonstrate a critical 
velocity for heating of the atom cloud |7| . Although the- 
oretical analysis of the GP equation has interpreted these 
results in terms of soliton and vortex production 0, Q , it 
is unclear whether quantum fluctuations arc important. 
These ideas lead us to pose quite fundamental and gen- 
eral questions about how the coherence of BECs degrades 
during supersonic transport past an obstacle, a defect in 
the potential, or more extended potentials, either peri- 
odic or disordered. 

In this paper, we address these questions by per- 
forming hundreds of three-dimensional numerical simu- 
lations of dilute gas BECs using the Truncated Wigner 
method [§, [l^ • This allows us to map out the flow of 
coherent and incoherent particles as we vary the BEC ve- 
locity and the nature of the obstacle or defect. We begin 
by studying the supersonic transport of a BEC around 
an impenetrable cylindrical obstacle. Previous analysis 
of the GP equation has shown that a Mach cone develops 
behind (or in the wake of) the obstacle, containing soli- 



tons and vortices 0, Q. In front of the obstacle, outside 
the Mach cone, a density modulation appears which is a 
form of Cherenkov radiation. Our calculations show that, 
in the presence of quantum fluctuations and strong inter- 
atomic interactions, this modulation exists only briefly, 
being almost completely destroyed at large times. Atoms 
are scattered incoherently into random directions, creat- 
ing quantum turbulence 0] and a crescent of incoherent 
atoms in front of the obstacle. 



These findings lead us to re-visit the classic laser 
stirring experiments 7], which showed heating above a 
critical velocity approximately 25% the bulk speed of 
sound. This effect has been attributed to vortex produc- 
tion 0, Q . Our simulations show that incoherent atoms 
may radiate from the oscillating laser but, in contrast to 
the vortex production, only at supersonic speeds since 
the Mach cone and Cherenkov radiation do not occur for 
subsonic transport. We quantify the incoherent scatter- 
ing by plotting the peak incoherent density as a function 
of BEC velocity, and hence demonstrate the threshold 
behaviour at the bulk speed of sound. 



Finally, we explore the supersonic transport of BECs 
through a disordered potential, which is lower in ampli- 
tude than the BEC chemical potential. Our results add 
to the interpretation of recent experiments which studied 
damped dipole oscillations in disordered potentials [ll|. 
Despite the BEC velocity reaching deep into the super- 
sonic regime in the experiments, hitherto no theoretical 
work has considered incoherent scattering. Theoretical 
work has been limited to studies of the GP equation, 
which has revealed coherent backscattering [l^ and soli- 
ton production , suggesting that there may be signif- 
icant heating and a reduction of the condensate fraction 
due to incoherent processes. Our simulations confirm 
this, revealing an accumulation of incoherent atoms as 
the condensate enters the disorder, and a corresponding 
suppression of the propagation of backscattered atoms. 



II. SYSTEM AND METHODOLOGY 

We consider a three-dimensional box, which has d 
mensions Lx, Ly and Lz, fiUed with a ^^Na BEC ; 
constant density no = 2.2 x 10^° m~'^, with per 
odic boundary conditions. We expand the wavefunctic 
^' {x, y, z, t), where t is time, over a basis of M plane- wa^ 
states, which are propagated in time using the fourtl 
order Runga-Kutta in the interaction picture algorith 
(RK4IP) [i, HI, [ll]. A flow is introduced by imposir 
a phase gradient in the x-direction. Due to the period 
boundary conditions, we are constrained by the conditic 
that the phase must change by 27rj over L^, where j is £ 
integer. The velocity is therefore = {h/m) x {2TTj /Lx). 
We investigate supersonic effects b y increasing Vx be- 
yond the bulk speed of sound Vg — \J h^noa/TrmP, where 
a — 2.9 nm is the s-wave scattering length and m is the 
mass of a single ^'^Na atom. 

In order to investigate supersonic flow of BECs un- 
der different conditions, we introduce barriers and de- 
fects in the potential. We minimise non-equilibrium ef- 
fects by varying the density across the barrier or de- 
fect such that the BEC chemical potential at t = is 
a constant throughout the grid. To avoid spurious ef- 
fects associated with the periodic boundary conditions 
our simulations never run for longer than a maximum 
time tmax = Lx/ (2vx)- 

Without quantum fluctuations, our scheme is equiva- 
lent to the projected GP equation [l4j. The GP equa- 
tion can describe the production of solitons, vortices, and 
Cherenkov radiation, but assumes complete coherence is 
always maintained. We include the effect of incoherent 
scattering in our model by adding random complex noise 
to each plane-wave mode. Heuristically, this simulates 
quantum vacuum fluctuations by adding classical fluctu- 
ations to the coherent field of the BEC's initial state. 
The amplitude of the noise has a Gaussian distribution, 
with an average value of half a particle per plane-wave 
mode, hence adding ~ M/2 "virtual particles" to the 
^ noLxLyLz "real particles" jl^. This technique is 
known as the Truncated Wigner method 0, ■ It mod- 
els the dynamics of incoherently scattered atoms, which 
are forbidden within the GP formalism. Specifically, 
the technique has been applied to colliding BECs [J], 
BEC refiection from a steep barrier [1], three-body re- 
combination processes [l3l, collapsing BECs atom 
interferometers [l^ , BEC diffraction 20] and optical lat- 
tices 21!]. 

By running multiple simulations with different initial 
random quantum fluctuations, we may calculate the dy- 
namics of coherent and incoherent atoms (2^ . Firstly, 
we compute the average density profile of the ensemble 
of simulations, then subtract the uniform density contri- 
bution M/ (2LxLyLz) from the virtual particles. This 
is the density profile of all atoms, coherent and incoher- 
ent. Secondly, we average ^' (x, y, z, t) over the ensemble 
of simulations, and compute the square modulus. This 
is the density profile of the coherent atoms. We may 
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FIG. 1: Density profiles in the z = plane (axes inset) of 
single simulations of a BEC flowing at Vx ~ l-5vs past an 
inpenetrable obstacle, whose position is indicated by a black 
cross, without (top row) and with (bottom row) quantum 
fluctuations, at t = [(a) and (e)], t = 0.40 ms [(b) and (f)], 
t = 0.64 ms [(c) and (g)], and t = 0.88 ms [(d) and (h)]. The 
field of view is 20/im x 20/im. Arrows in (d) and (h) indicate 
direction of quantized circulation around vortex cores. 



then calculate the profile of incoherent atoms by simply 
subtracting the coherent density profile from the total 
density profile. 



III. FLOW PAST AN INPENETRABLE 
OBSTACLE 

A. Single trajectory results 

The top row of Fig. [1] shows the well-known mean-field 
picture of Cherenkov radiation and vortex production as 
a BEC fiows past an inpenetrable obstacle at a supersonic 
speed. These images are density profiles in the z — 
plane, calculated using the bare GP equation without 
quantum fiuctuations. Hence, at i = 0, the BEC has a 
uniform density away from the obstacle [Fig. (TJa)]. The 
obstacle, whose position is indicated by the black cross in 
Fig. [IJa), is a cylinder of diameter 1 micron and length 
Lz, with its axis running parallel to the z-axis. Within 
the obstacle the atom density is set to zero. The imposed 
condensate velocity Vx = 1.4'i;s- 

The Cherenkov radiation is observed as density mod- 
ulations ahead of (to the left of) and curving around the 
obstacle [Fig. [Tl]b)-(d)], forming a cone shape known as 
the Mach cone. The angle 6 between the flow direction 
and the Mach cone line is given by Q 



(Vs/Vx) 



(1) 



As t increases from 0.40 [Fig.IUJb)] to 0.88 ms [Fig.fHd)], 
the radiation gradually extends away from the obstacle. 
Vortices are seeded in the wake of (to the right of) the 
obstacle within the Mach cone [Fig. [TJc) and (d)] , and 
subsequently travel away in the direction of condensate 
flow. Two vortices are encircled by arrows in Fig. md), 
indicating the direction of circulation. 



The bottom row of Fig. [T] shows the correspondh 
images including quantum fluctuations [Fig. [He)]. I 
t = 0.40 ms [Fig. [ijf)], the Cherenkov radiation can I 
observed as it begins to form. However, it is not visib 
for t > 0.50 ms [Fig.[T]Jg) and (h)], and has been replace 
by a turbulent flow, including vortices [encircled by tl 
white arrows in Fig. [ijh)] which are completely absc] 
from the corresponding GP picture [Fig. [T]^d)]. This tu 
bulent flow has been referred to as "quantum turbulence 
in the context of colliding BECs [1] . The vortices forme 
in the wake of (to the right of) the obstruction are una 
fected by quantum fluctuations, and consequently appe; 
in identical locations to those in the GP simulations. T'w 
examples are encircled by the black arrows in Fig. [Hhj, 
indicating the direction of circulation. 



B. Trajectory ensemble results 

We may gain more insight by averaging multiple sim- 
ulations with different initial quantum fluctuations. We 
calculate the average total density of one hundred simu- 



lations within a cube of dimensions Lr 



20 



/im, then integrate along the z-axis to obtain a two- 
dimensional density profile in the x-y plane, equivalent to 
an experimental absorption image. The result is shown in 
the top row of Fig. O Since the quantum fluctuations are 
random, the initial profile of the averaged atom density 
appears uniform [Fig. El^a)], except within the obstacle 
(whose position is indicated by the black cross) where 
the density is zero. These averaged images confirm that 
the Cherenkov radiation visible a.tt = 0.40 ms [Fig.[2lb)] 
is strongly supressed hy t > 0.50 ms [Fig. [2lc) and (d)], 
in contrast to the clear density modulations found in the 
absence of quantum fluctuations [Fig. [IJc) and (d)]. As 
noted above, the production of vortices within the Mach 
cone [Fig.[2lb)-(d)] is identical to that in the simulations 
excluding quantum fluctuations [Fig. [T](b)-(d)]. 

The corresponding incoherent atom density profiles are 
shown in the lower row of Fig. [51 At t ~ the simu- 
lations are completely coherent, and hence the incoher- 
ent atom density profile is zero [Fig. [2l[e)]. At t = 0.40 
ms, when the Cherenkov radiation is still clearly visi- 
ble [Fig. [D(b)], the incoherent density remains very low 
[Fig. [2[f)]. However, as the Cherenkov radiation begins 
to fade [Fig. [ll^c)], a crescent-shaped concentration of in- 
coherent atoms appears ahead of (to the left of) the ob- 
struction [Fig. [2tg)]. The shape of the incoherent cloud 
reflects the regions of intense Cherenkov radiation. As 
t increases to 0.88 ms, the crescent of incoherent atoms 
gradually extends around the Mach cone and very slowly 
away from the cone [Fig. [H^h)], much slower than the 
Cherenkov radiation. These results show that the forma- 
tion of Cherenkov radiation is suppressed by incoherent 
scattering in dense BECs with strong inter-atomic inter- 
actions, causing heating and a reduction in the conden- 
sate fraction. 

Cherenkov radiation can be distinguished from inco- 
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FIG. 2: Simulated absorption images in x — y plane (axes in- 
set) averaged over one hundred simulations of a BEC moving 
past an inpenetrable obstacle, whose position is indicated by 
a black cross, for Vx ~ l-5vs, at t = [(a) and (e)], t = 0.40 ms 
[(b) and (f)], t = 0.64 ms [(c) and (g)], and t = 0.88 ms [(d) 
and (h)]. Top (bottom) row shows total (incoherent) density. 
The field of view is 20/xm x 20^m. Arrows in (d) indicate 
direction of quantized circulation around vortex cores. 



herently scattered atoms by studying their dynamics. As 
noted above, the center-of-mass of the incoherent atoms 
is approximately stationary (see Fig. [2]). This is because 
these atoms have been scattered into random directions 
in the regions of intense Cherenkov radiation, and hence 
some atoms are travelling directly towards, or directly 
away from, the observer. In addition, multiple scattering 
events may occur, hence suppressing the propagation of 
the incoherent atoms. In contrast, the Cherenkov radia- 
tion is coherently emitted from the obstruction, travelling 
in a direction perpendicular to the surface of the Mach 
cone, as shown in Fig. [1] As an example, recent experi- 
ments [l| have reported a shock-wave which propagates 
away from the obstacle during supersonic transport. This 
propagation identifies the shock-wave as Cherenkov radi- 
ation; incoherently scattered atoms would remain close 
to the obstacle. 

In Fig. [31 we quantify the incoherent scattering by plot- 
ting the averaged peak incoherent atom density TL^ as a 
function of Vx t = 0.4tinax- The curve shows that rii 
rises abruptly from zero when Vx exceeds Vs , approaching 
a maximum value, set by no, for Vx ~ Sw^. This threshold 
behavior occurs because the Mach cone and Cherenkov 
radiation do not form for flow speeds below the speed of 
sound ^ {9 has a maximum value of 7r/2 for Vx = Vs, see 
Eq. [1]). In contrast, vortex shedding can occur for flow 
speeds below the speed of sound, as has been observed 
experimentally 0. To illustrate this, the inset in Fig. [31 
shows the density proflle in the 2 = plane at t — 1. 1 ms 
(= O.Timax), for Vx — O.Svg, without quantum fluctua- 
tions. No Mach cone or Cherenkov radiation is present, 
but two vortex anti-vortex pairs (encircled by arrows to 
indicate the direction of circulation) have been shed from 
the obstacle. Vortices can be generated for flow speeds 
below the speed of sound because the relative velocity 
between the obstacle and the BEC may exceed the speed 
of sound locally Q- 




FIG. 3: Hi versus Vx &t t = 0.4tmax for a BEG flowing past 
an inpenetrable obstacle. Inset shows the density proflle in 
the z = plane (axes inset) for Vx = 0.8vs, without quantum 
fluctuations. The position of the obstacle is indicated by the 
black cross. The fleld of view is lOfim x 12^m. Arrows in 
inset indicate direction of quantized circulation around vortex 
cores. 
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FIG. 4: Simulated absorption images in x ~ y plane (axes 
inset) averaged over one hundred simulations of a BEG stirred 
by laser, whose position is indicated by a black cross, with 
vo = 1.8«s, at t = [(a) and (e)], t = 0.36 ms [(b) and (f)], 
t = 0.58 ms [(c) and (g)], and t = 0.80 ms [(d) and (h)]. Top 
(bottom) row shows total (incoherent) density. The field of 
view is 20/im x 20/im. Arrows in (c) indicate direction of 
quantized circulation around vortex cores. 



IV. LASER STIRRING 

The findings in the previous section motivated us to re- 
visit the classic laser stirring experiment, in which heat- 
ing was observed above a critical velocity . Theoretical 
work has attributed this heating to soliton and vortex 
production 0, , but incoherent scattering may play a 
significant role. We therefore replace the stationary ob- 
stacle in the previous section with an obstacle oscillating 
in the x-direction with speed Vo = ±1.8us and amplitude 
A = 4.0 /zm, and hence period to — 4:A/vo — 1.2 ms. 
The imposed phase gradient is also removed so the BEC 
is initially stationary. 

The top row of Fig. 0] shows the simulated absorption 
images generated from one hundred laser drag simula- 
tions. The BEC density is uniform at i = 0, except near 
the black cross which indicates the position of the laser 
[Fig.lU^a)]. Initially, the BEC is completely coherent, and 
hence the incoherent density is zero [Fig. dfe)]. As the 
laser moves from right to left [Fig. |3]Ja)-(c)], Cherenkov 
radiation moves ahead of (to the left of) the laser, and 
vortices are generated in its wake (to its right). A cres- 
cent of incoherent atoms develops [Fig. HJ^f) and (g)], 
moving from right to left ahead of the obstacle. We can 
again interpret this as incoherent scattering due to the 
Cherenkov radiation. 

For t < to/ 2 — 0.58 ms, the dynamics are identical 
to those produced by the static obstacle in the previous 
section: we have simply changed the observer's frame- 
of-reference. But after the laser reverses its direction 
of motion, the two systems are no longer equivalent. 
The Cherenkov radiation and the crescent of incoherent 
atoms continue to propagate in the original rest frame 
of the laser, whilst the laser travels in the opposite di- 
rection [Fig.lD^d) and (h)]. The supersonic motion of the 
Cherenkov radiation and incoherent atoms causes further 
incoherent scattering, broadening the incoherent region 



and raising the peak incoherent density. Meanwhile, the 
laser moves in the positive z-drection, and we now begin 
to observe the formation of a second crescent of incoher- 
ent atoms. 

Incoherent scattering and subsequent quantum turbu- 
lence will contribute to the heating observed experimen- 
tally when a laser oscillates in a BEC Q. However, as 
shown in Fig. [31 this additional heating occurs only when 
Vo > Vs, becoming significant for Vo > l-5vs- The experi- 
ments focused instead on the threshold for heating, which 
occurs at a critical velocity of approximately 0.25us. In 
this threshold regime, the heating will be solely due to 
coherent vortex production [2|, Q , because Cherenkov ra- 
diation does not occur below the bulk speed of sound. 
Consquently, these findings do not significantly change 
the interpretation of the experiments [71 . 



V. DISORDERED POTENTIALS 

The previous sections have discussed incoherent scat- 
tering from inpenetrable barriers, within which the po- 
tential exceeds the BEC chemical potential. However, 
significant scattering is also possible from small poten- 
tial defects, lower in amplitude than the BEC chemical 
potential. To investigate this problem we generate a ran- 
dom potential with a peak amplitude of half the chemi- 
cal potential. We set the potential to zero in the region 
X < in order to highlight the effect of backscattering 
from the random potential. The initial density of the 
BEC is adjusted so that its chemical potential is a con- 
stant. Consequently the simulated absorption image at 
t = is uniform for x < [Fig. EJa), x = is indicated 
by the inset axes], and varies smoothly over the random 
potential for a; > 0. 

The top row in Fig. [5] shows the dynamics of the 
BEC in the disordered potential without quantum 
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FIG. 5: Top row: simulated absorption images in x — y plane 
(axes inset), without quantum fluctuations, of a BEG in a 
disordered potential, for Vx = l.Sns, at t — (a), t = 1.0 ms 
(b), t = 2.0 ms (c), and t = 3.0 ms (d). Middle (bottom) 
row: corresponding images including quantum fluctuations, 
averaged over one hundred simulations, showing total (inco- 
herent) atom density. The field of view is 20/im x 20/im. 



fluctuations. The random potential causes a back- 
scattering [l^, which can be observed as the interfer- 
ence pattern close to the boundary between the flat and 
random potential regions at cc = [Fig. Eljb)]. As time 
increases, this interference pattern increases in amplitude 
and extends away from the boundary in the negative a;- 
direction [Fig.l^l^c) and (d)]. 

If quantum fluctuations are included, the initial density 
profile is identical [Fig. [Sl^e)], and the corresponding in- 
coherent density is zero [Fig. . Backscattering begins 
to occur [Fig. O^f)] as before, but, rather than increas- 
ing in amplitude and extending away from the boundary 
region, the interference pattern begins to fade, and re- 
mains close to X = [Fig. E^g) and (h)] . This is due to 
incoherent scattering between the coherently backscat- 
tered atoms and the incoming atoms approaching the 



disordered potential. Since the atoms are scattered into 
random directions, they remain close to the boundary re- 
gion [Fig. [5]^j)-(l)], rather than propagate away like the 
coherent backscattered atoms. This process introduces 
an additional contribution, which is not described within 
the GP equation, to the damping of BEG motion ob- 
served experimentally TH . There are also several regions 
of high incoherence within the random potential, due to 
backscattering from local peaks in the potential. 



VI. CONCLUSIONS 

In summary, we have shown that supersonic transport 
of BEGs can generate significant incoherent scattering. 
This causes a reduction of the condensed fraction, and 
hence an additional contribution to the loss of superflu- 
idity and the supression of BEG motion, which is typi- 
cally observed in experiment as heating and a damping of 
dipole oscillations 0, [U . We have mapped out the dy- 
namics of the incoherent atoms by performing hundreds 
of three-dimensional simulations, and related our results 
to relevant experiments. 

In principle, our predictions can be tested experimen- 
tally, but in practice it may be difficult to separate heat- 
ing due to coherent soliton or vortex production from 
heating due to incoherent scattering. It may be possible 
to exploit the drop in the condensed fraction as a result of 
incoherent scattering by, for example, applying a Bragg 
pulse to create fringes in the cloud. Due to loss of phase 
coherence, there will be a reduction in the fringe visibil- 
ity [11]. In contrast, coherent backscattering and vortex 
production will generate a reproducable fringe pattern. 
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